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Abstract 

Background: The release of radioactive materials due to the Fukushima nuclear accident has raised concern 
regarding the biological impacts of ingesting radioactively contaminated diets on organisms. We previously 
performed an internal exposure experiment in which contaminated leaves collected from polluted areas were 
fed to larvae of the pale grass blue butterfly, Zizeeria maha, from Okinawa, which is one of the least polluted 
localities in Japan. Using the same experimental system, in the present study, we further examined the effects 
of low-level-contaminated diets on this butterfly. Leaves were collected from two localities in Tohoku (Motomiya 
(161 Bq/kg) and Koriyama (1 1 7 Bq/kg)); two in Kanto (Kashiwa (47.6 Bq/kg) and Musashino (6.4 Bq/kg)); one in Tokai 
(Atami (2.5 Bq/kg)); and from Okinawa (0.2 Bq/kg). In addition to the effects on the first generation, we examined the 
possible transgenerational effects of the diets on the next generation. 

Results: In the first generation, the Tohoku groups showed higher rates of mortality and abnormalities and a smaller 
forewing size than the Okinawa group. The mortality rates were largely dependent on the ingested dose of caesium. 
The survival rates of the Kanto-Tokai groups were greater than 80%, but the rates in the Tohoku groups were much 
lower. In the next generation, the survival rates in the Tohoku groups were below 20%, whereas those of the Okinawa 
groups were above 70%. The survival rates in the second generation were independent of the locality of the leaves 
ingested by the first generation, indicating that the diet in the second generation was the determinant of their survival. 
Moreover, a smaller forewing size was observed in the Tohoku groups in the second generation. However, the 
forewing size was inversely correlated with the cumulative caesium dose ingested throughout the first and 
second generations, indicating that the diet in the first generation also influenced the forewing size of the 
second generation. 

Conclusions: Biological effects are detectable under a low ingested dose of radioactivity from a contaminated diet. 
The effects are transgenerational but can be overcome by ingesting a non-contaminated diet, suggesting that at least 
some of the observed effects are attributable to non-genetic physiological changes. 

Keywords: Fukushima nuclear accident, Ingestion, Internal exposure, Low-dose exposure, Pale grass blue butterfly, 
Radioactive contamination, Transgenerational effect 



Background 

The collapse of the Fukushima Dai-ichi Nuclear Power 
Plant (NPP) due to the Great East Japan Earthquake on 11 
March 2011 resulted in the release of a massive amount of 
radioactive materials into the surrounding environment. 
This large-scale environmental pollution has impacted the 
lives of organisms living in the polluted areas. A field 
study revealed that the number of butterflies decreased in 
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the polluted areas [1,2], and morphological abnormalities 
have been detected in aphids from Fukushima [3]. How- 
ever, such studies examining the biological impacts of the 
accident are still scarce. 

Our group has been studying the development and 
evolution of butterflies in response to environmental 
stress [4]. One of our favourite butterflies is a small ly- 
caenid butterfly, the pale grass blue, Zizeeria maha 
(Lepidoptera, Lycaenidae), which is found throughout 
Japan, except in Hokkaido [4-6]. Using this butterfly as 
an experimental system, we previously reported various 
morphological abnormalities and a smaller forewing size 
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in adults collected from highly polluted areas [7-9]. We 
showed that some abnormal traits, including aberrations 
in the wing colour pattern, were heritable, suggesting 
genetic damage introduced by the accident [7-9]. We be- 
lieve that the proposed genetic damage was introduced 
mainly by the initial exposure of the butterflies immedi- 
ately after the explosion [9]. We reproduced these field- 
based results through an external exposure experiment in 
which butterfly larvae and pupae were raised near an arti- 
ficial caesium radiation source as well as through an in- 
ternal exposure experiment in which larvae consumed 
contaminated host plant leaves collected from polluted 
areas [7-9]. High mortality and abnormality rates were de- 
tected [7], and a small forewing size (i.e., body size) was 
also observed [7], similar to the outcome of the food stress 
in butterflies [10]. 

In a previous study, we collected host plant leaves of 
Oxalis corniculata (in the summer of 2011) from four 
relatively highly polluted localities in the Tohoku district 
(Iitate-montane, Iitate-flatland, Fukushima, and Hirono), 
in addition to the control locality, Ube (Figure 1) [7]. We 
further quantified caesium radioactivity in pupae that 
ingested contaminated leaves and determined the mathem- 
atical relationships between the ingested dose of caesium 
and the resultant rates of mortality and abnormalities [11]. 



We found that the mortality and abnormality rates in- 
creased sharply under low doses of caesium ingestion, 
following a power function. However, as the leaves show- 
ing the lowest concentration recorded in that study were 
collected from Hirono (1,452 Bq/kg), leaves from other lo- 
calities with much lower contamination levels should be 
examined to better understand the effects of low-level ra- 
diation ingestion. 

In the present study, we collected leaves from host 
plants of the pale grass blue butterfly from six localities 
showing relatively low-level contamination: Motomiya, 
Koriyama (both in Fukushima Prefecture in the Tohoku 
district; collected in the fall of 2012), Kashiwa and 
Musashino, (in Chiba Prefecture and Tokyo, respectively, 
in the Kanto district; collected in the summer of 2012), 
Atami (in Shizuoka Prefecture in the Tokai district; col- 
lected in the summer of 2012) and Okinawa (one of the 
farthest districts from the Fukushima Dai-ichi NPP in 
Japan; collected in the summer and fall of 2012) (Figure 1; 
Table 1). These host plant leaves were given to larvae from 
Okinawa, which were not affected by the nuclear accident 
to any detectable degree, under our standard rearing 
conditions [12]. We evaluated the biological impacts of 
various levels of contamination on this butterfly in this 
manner. Furthermore, we performed the internal 




.itate Village 
1 •Ji NPP 

ono Town 



Musashino City • ^KasTiiwa City 
(/Atami City 



Nishihara Town (Okinawa) 



200 km 



Figure 1 Collection localities for host plant leaves. The localities examined in the present study are shown in blue letters and those examined 
in the previous studies [7-10] in black letters. The Fukushima Dai-ichi Nuclear Power Plant is indicated as NPP in red. 
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Table 1 Information for the collection localities for host plant leaves 





Collection site 


Collection date 


Distance from the 
NPP [km] (average) 


Ground radiation 
[uSv/h] (average) 


Okinawa 


1 Senbaru, Nishihara, Okinawa Prefecture, University of the Ryukyus 


20, 30 July 2012 
2-22 October 2012 


1760 


0.05 (8 April 2012) 


Atami 


1 164 Izusan, Atami City, Shizuoka Prefecture 


25 July 2012 
4 August 2012 


311 


0.07 


Musashino 


2 Gotenyama, Musashino City, Tokyo 


24 July 2012 
3 August 2012 
20 October 201 2* 1 


232 


0.12 


Kashiwa 


6 Kashiwanoha, Kashiwa City, Chiba Prefecture 


23 July 2012 
7 August 2012 
5 October 201 2* 1 


196 


0.47 


Koriyama 


Tomita, Koriyama City, Fukushima Prefecture 


4, 14 October 2012 
1 November 2012 


61 


1.13 


Motomiya 


Arai, Motomiya City, Fukushima Prefecture 


4, 14 October 2012 
1 November 2012 


59 


1.42 



* 1 Leaves were collected only for activity measurements. 

exposure experiment not only in the first (F^ generation, 
but also in the second (F 2 ) generation to examine possible 
transgenerational effects. 

Results 

Effects on the Ft generation 

Previous studies showed that an average larva has a 
body weight of 0.035 g and consumes 0.388 g of leaves 
throughout its life on average [11]. We obtained the con- 
centrations of caesium activity in leaves collected from 
the six indicated localities (Table 2), and we calcu- 
lated the amount of caesium activity ingested by each 
feeding group of larvae, together with the resultant 
mortality and abnormality rates (Table 3). The mor- 
tality rate and the abnormality rate were defined as per- 
cent dead or percent abnormal and dead (and hence, total 
abnormality rate), respectively. We then examined the 
dose-response relationship between the ingested caesium 

Table 2 Caesium radioactivity in host plant leaves* 1 



dose and the mortality rate (Figure 2a). The mortality 
rate (y) increased linearly in accordance with an in- 
crease of the caesium dose {x) with the regression 
model, y = 0.60 (±0.21) x + 10.21 (±6.23) {R 2 = 0.659, df= 
5, F = 7.71, p = 0.0497). Remarkably, the mortality rate of 
the Koriyama group (40.9 mBq per larva) was 53%. 
The dose-response relationship for abnormality rates was 
almost identical to that for mortality rates (not shown). 

Because surviving adults are not necessarily healthy re- 
garding morphological traits, we calculated relative nor- 
mality values for adults in reference to the control 
Okinawa groups as an indication of the robustness of 
each feeding group (Figure 2b). To this end, the ratio of 
the number of morphologically normal (non-abnormal) 
adults to the original number of larvae was calculated, 
and the ratios in the control Okinawa groups were set at 
100%. As expected, the robustness distribution (y) ap- 
peared to show a non-significant inverse response 



Activity of 137 Cs (±SE) [Bq/kg] 



Activity of 134 Cs (±SE) [Bq/kg] 



Total for 137 Cs and 134 Cs [Bq/kg] 



Okinawa (for 


0.1 2 ±0.02 


0.06 ± 0.02 


0.1 8 ±0.03 


Atami (for F,) 


1 .48 ± 0.06 


1 .05 ± 0.06 


2.53 ± 0.08 


Musashino (for F,) 


3.82 ± 0.07 


2.56 ±0.07 


6.38 ±0.10 


Kashiwa (for F,) 


28.40 ± 0.33 


19.1 7 ±0.32 


47.57 ±0.46 


Koriyama (for F,) 


71.87 ±0.40 


45.34 ±0.37 


117.21 ±0.54 


Motomiya (for F,) 


98.1 9 ±0.54 


62.41 ±0.49 


160.60 ±0.73 


Okinawa (for F 2 ) 


0.11 ±0.02 


0.06 ±0.02 


0.1 8 ±0.03 


Koriyama (for F 2 ) 


71.73 ±0.40 


44.07 ± 0.36 


11 5.80 ±0.54 


Motomiya (for F 2 ) 


98.00 ± 0.53 


60.65 ± 0.48 


158.65 ±0.72 



#1 Activity levels were calculated assuming that the larvae ate all of the leaves required for their subsequent growth at once on the first day of ingestion of the 
contaminated leaves and that 1 37 Cs and 1 34 Cs were released at a 1:1 activity ratio on 1 5 March 201 1 in a single burst from the Fukushima Dai-ichi NPP. 
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Table 3 Ingested caesium dose and mortality and abnormality rates in the F n 


Ud Iwl CI llwl 1 




The f^ leaves from: 


Number of larvae (n) 


Ingested caesium dose [mBq] 


Mortality rate [%] 


Total abnormality rate [%] 


Okinawa 


308 (Summer = 67, Fall = 241) 


0.0629 + 0.0011 


8.05* 1 


8.25* 1 


Atami 


55 


0.883 ± 0.003 


9.10 


9.10 


Musashino 


54 


2.23 ±0.004 


13.0 


13.0 


Kashiwa 


55 


16.6 ±0.02 


16.4 


16.4 


Koriyama 


249 


40.9 ± 0.02 


53.0 


54.2 


Motomiya 


256 


56.1 ±0.03 


31.2 


32.0 



* 1 Mean values for the summer and fall results. 



relationship with the caesium dose (x) with the regres- 
sion model, y = -0.53 (±0.24) x + 89.07 (±6.23) (R 2 = 
0.555, df=5, F=4.99, p = 0.089). When the Motomiya 
sample was excluded as an outlier, the regression 
equation becomes y = -1.03 (±0.17) x + 92.37 (±3.41) 
{R 2 = 0.923, df= 4, F= 35.95, p = 0.0093), showing a highly 
significant dose-response relationship and being consist- 
ent with the results found for the mortality rates. 



We measured the forewing size of the surviving adults 
from each feeding group, and these data were plotted 
against the ingested caesium dose (Figure 3). The fore- 
wing size fluctuated at low doses, but it decreased at 
relatively high doses in both sexes. Statistically signifi- 
cant differences were obtained between the Okinawa 
and Koriyama groups both in males (t = 2.22, df= 137, 
p = 0.028, Welch's f-test) and in females (£=3.17, df = 98, 
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Figure 2 Effects of the ingested caesium dose on the F, generation, (a) Mortality rate, (b) Relative normality rate (robustness). 
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p = 0.0020, Welch's i-test) and between the Okinawa and 
Motomiya groups both in males (t = 3.31, df= 166, p = 
0.0011, Student's f-test) and in females (f=3.51, df= 155, 
p = 0.0006, Welch's f-test). A linear regression model 
between the ingested caesium dose (x) and the fore- 
wing size (y) did not show statistically significant rela- 
tionship in males with the equation, y = -0.00083 
(±0.00034) *+ 1.0215 (±0.0099) (R 2 = 0.598, df=5, F = 
5.95, p = 0.071), but it showed a significant relation- 
ship in females with the equation, y = -0.00076 
(±0.00017) x+ 1.0098 (±0.0050) (R 2 = 0.827, df=5, F = 
19.14, p = 0.012). This sexual difference may originate 
from the simple fact that females have larger wings in size. 

We generated developmental survival curves (which in- 
dicate changes of survival rate, i.e., percent survival, over 
developmental stages) for the six feeding groups (Figure 4). 
These survival curves showed stage-dependent changes of 
the percentage of the surviving individuals. The survival 
curves for the five groups examined in a previous study 



[7] were also included in the figure for comparison. Most 
of the curves were largely parallel with one another, and 
as a whole, these survival curves were statistically different 
from one another {df= 11, % 2 = 300, p< 0.0001, log rank 
test; df= 11, x 2 = 292, p < 0.0001, Wilcoxon rank test). 

We noted that the Motomiya and Koriyama curves 
showed a great deal of mortality in the larva-prepupa 
period. The curves were segregated into three clusters of 
groups: a high survival cluster (Okinawa, Ube, Atami, 
Musashino, and Kashiwa), an intermediate survival clus- 
ter (Hirono and Motomiya), and a low survival cluster 
(Iitate-montane, Koriyama, Fukushima, and Iitate- 
flatland). These three clusters were statistically different 
from one another (df= 3, % 2 = 286, p < 0.0001, log rank 
test; df=3, / 2 = 276, p< 0.0001, Wilcoxon rank test). 
However, the concentrations of caesium activity re- 
corded in the leaves showed great variation within a 
given cluster. This intra-cluster variation most likely oc- 
curred because these groups were implemented using 
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Figure 3 Relationship between the ingested caesium dose and wing size in the F n generation. The relative mean wing size is shown, 
together with standard error bars. The size in the Okinawa group was set as 1.00. (a) Male forewings. The number of individuals that were 
subjected to measurement is as follows: Okinawa (n = 97), Atami (n = 25), Musashino (n = 1 9), Kashiwa (n = 25), Koriyama (n = 84), and Motomiya 
(n = 71). (b) Female forewings. The number of individuals that were subjected to measurement is as follows: Okinawa (n = 99), Atami (n = 23), 
Musashino (n = 27), Kashiwa (n = 21), Koriyama (n = 63), and Motomiya (n = 87). 
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three different genetic lines of this butterfly species in 
different time periods, suggesting intraspecific variation 
regarding sensitivity to a radioactively contaminated diet. 
Accordingly, a set of groups that were assessed within 
the same time period were found to be comparable to 
one another as they come from the same genetic line 
(see Methods). 

We discovered various morphological abnormalities in 
the surviving adults (Figure 5). The severe and rare ab- 
normalities shown in Figure 5 might imply the effects of 
a contaminated diet. Only three individuals in the group 
that consumed the control Okinawa leaves showed very 
minor morphological abnormalities, i.e., folded wings 
and a deformed leg, which were not comparably severe 
with the abnormal individuals from the other F x groups. 

Effects on the F 2 generation 

The F 2 generation was obtained from the morphologic- 
ally healthy (non-abnormal) Y 1 adults that consumed 
leaves from Okinawa, Koriyama, or Motomiya. Based on 
the detected caesium activity in the leaves, the ingested 
caesium doses were calculated (Table 4), together with 
the resultant mortality and abnormality rates (Table 5), 
and the dose-response relationship between the ingested 
dose and the mortality rate was examined. The mortality 
rate (y) appeared to be dependent on the caesium dose 
(x) ingested by the F 2 generation with the regression 
model, y = 1.445 (±0.088) *+ 17.04 (±3.24) (R 2 = 0.982, 
df= 6, F = 267, p < 0.0001), but not significantly by the Fj 



generation (R 2 = 0.00675, df=6, F= 0.034, p = 0.86) 
(Figure 6a, b). 

Next, we produced developmental survival curves for 
the F 2 generation (Figure 7). As a whole, these survival 
curves were highly different from one another (df= 6, 
X 2 = 818, p < 0.0001, log rank test; df= 6, / 2 = 824, p < 
0.0001, Wilcoxon rank test). The curves were clearly seg- 
regated into two clusters: high and very low survival clus- 
ters. These two groups were statistically very different 
from each other (df= 1, / 2 = 758, p < 0.0001, log rank test; 
df= 1, x 2 = 774, p < 0.0001, Wilcoxon rank test). Again, it 
was apparent that the determinant of this segregation was 
the F 2 leaves and not the Fi leaves. The remarkable segre- 
gation into the two clusters indicated the importance of 
the F 2 diet. 

We compared the forewing sizes of the surviving adults 
(Figure 8a,b). The group that consumed the Koriyama 
leaves in both the ¥ 1 and F 2 generations showed very small 
forewings in both sexes compared to the group that con- 
sumed the Okinawa leaves in both generations (£ = 8.44, 
df= 74, p < 0.0001, Student's i-test, for males; t = 5.03, df 
= 13, p = 0.0002, Welch's f-test, for females). Interestingly, 
the cumulative ingested caesium dose throughout the Fi 
and F 2 generations showed a higher correlation with the 
forewing size than the ingested dose in the F 2 generation 
alone in both sexes, but especially in females (Table 6; 
Figure 8c). The female regression line between the cumu- 
lative ingested dose (x) throughout the Fj and F 2 genera- 
tions and the forewing size (y) was obtained as 
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Figure 4 Survival curves for the F, generation. The six groups examined in the present study and the five groups examined in previous 
studies [7-10] are plotted together. The survival curves were roughly segregated into three clusters of the experimental groups. The concentrations of 
caesium activity in the ingested host plant leaves are indicated, together with the time period of the experiments. Groups that were examined in the 
summer of 2012 are indicated by brownish colours with circular symbols. Groups examined in the fall of 2012 are indicated by purplish colours with 
square symbols. Groups examined in the summer of 201 1 in previous studies [7-10] are indicated by greenish colours with triangular symbols. 
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y = -1.70 (±0.31) * + 1.259 (±0.020) (R 2 = 0.857, df=6, 
F= 29.99,^ = 0.0028). 

As observed in the Fi generation, various morpho- 
logical abnormalities were detected in the surviving F 2 
adults (Figure 9). Considering the very high mortality 
and abnormality rates recorded, it is likely that many of 
the recorded effects are attributable to the contaminated 
diets ingested by the butterflies. 



Table 4 Ingested caesium dose in the F 2 generation 



F 2 leaves from: 


Number of 


larvae (n) 


Ingested caesium 










dose [mBq] 


Okinawa 


522 


(=170 + 


1 28 + 224) 


0.0629 + 0.0011 


Koriyama 


290 


(=158 + 


132) 


40.4 ± 0.02 


Motomiya 


516 


(=158 + 


358) 


55.4 ± 0.03 



Discussion 

In the present study, we performed an internal exposure 
experiment using contaminated host plant leaves. The 
host plant leaves were collected from the localities show- 
ing low-level contamination, including two localities in 
Tohoku (Koriyama and Motomiya), two in Kanto (Kashiwa 
and Musashino), and one in Tokai (Atami). A similar ex- 
periment was performed previously [7], and the present 
study therefore serves to reproduce and validate the previ- 
ous results. Furthermore, the present study extends the 
previous study by using the leaves that were less contami- 
nated than those used in the previous study. 

In the present study, low-dose effects were clearly de- 
tected in the butterflies in the form of high mortality 
and abnormality rates, low normality rates (robustness), 
low survival rates, and a small forewing size in the ¥ 1 
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Table 5 Mortality and abnormality rates in the F 2 
generation* 1 



F 2 leaves 
F^ leaves 


Okinawa (F 2 ) 


Koriyama (F 2 ) 


Motomiya (F 2 ) 


Okinawa (FO 


MR 


= 16.5%; 


MR-- 


= 79.1%; 


MR- 


= 89.9%; 




AR-- 


= 1 8.8% 


AR- 


= 79.7% 


AR = 


= 90.5% 




{n = 


170) 


{n = 


158) 


(n = 


158) 


Koriyama (Fi) 


MR- 


= 10.2%; 


MR'- 


= 78.8%; 


Not 


tested. 




AR = 


= 1 5.6% 


AR- 


= 87.9% 








{n = 


128) 


{n = 


132) 






Motomiya (F^ 


MR- 


= 22.8%; 


Not 


tested. 


MR-- 


= 99.2%; 




AR-- 


= 29.9% 






AR = 


= 99.2% 




(n = 


224) 






(n = 


358) 



mortality rate; AR, total abnormality rate; n, number of larvae examined. 



generation of the Tohoku groups. It was somewhat sur- 
prising to find that the Koriyama group, which ingested 
40.9 mBq of radioactive caesium during the larval stage 
from contaminated-contaminated leaves showing an ac- 
tivity of 117 Bq/kg, displayed mortality and abnormality 
rates greater than 50%, although the Motomiya group 
exhibited somewhat lower rates. These results suggest 
that low-dose ingestion of approximately 100 Bq/kg may 
be seriously toxic to certain organisms. Even among the 
surviving adults, a small forewing size was detected in 
the Koriyama and Motomiya groups, confirming the re- 
sults concerning mortality rates, abnormality rates, nor- 
mality rates, and survival rates. 

For the Kanto and Tokai localities, which presented rela- 
tively low contamination levels ranging from 2.53 Bq/kg to 
47.57 Bq/kg, both the mortality (survival) rates and the ab- 
normality (normality) rates appeared to be dose dependent. 
These Kanto-Tokai groups were tested using the same 
genetic lines at the same time, which were comparable to 
one another. Therefore, the contaminated leaves from the 
Kanto-Tokai districts may have had small but detectable 
effects on the fitness of the butterflies. On the other hand, 
the forewing size of males in the Kanto-Tokai groups ap- 
peared to be larger than that of the Okinawa group, despite 
the possible adverse effects detected for other measures 
(Figure 3). The reasons for this size increase are not clear. 
The biological responses to very low doses may simply be 
different from the responses to relatively high doses. 

Interestingly, the survival curves obtained for the dif- 
ferent localities were segregated into three clusters: high, 
intermediate, and low survival clusters, not only among 
the groups from the previous study [7], but also in the 
groups of the present study, despite the large difference 
in caesium activity levels between the previous and 
present studies. In the previous study [7], the half mor- 
tality rate was calculated as 1.9 Bq per larva, corresponding 
to constant leaf consumption of 4.9 x 10 3 Bq/kg 



throughout the larval stage. These previous values were 
more than 40 times as high as the values we obtained for 
Koriyama in the present study. As a result, within a given 
cluster of survival curves, the levels of leaf contamination 
were highly variable (see Figure 4). However, within a 
given set of groups that were treated within the same time 
period (hence, coming from the same genetic line), the 
survival rates were largely dependent on the caesium ac- 
tivity in the leaves. These results may originate simply 
from a dependence of radiation sensitivity on the genetic 
background of butterfly lines. That is, direct comparison 
between different butterfly lines may be difficult, even if 
the experimental set-ups are virtually identical with the 
exception of their genetic makeup. Therefore, our results 
suggest the existence of intraspecific individual variation 
of sensitivity to low-level radiation, which is consistent 
with the existing literature [13]. 

The biological effects of ingesting the contaminated 
diets were more severe in the F 2 generation, despite the 
normal morphological appearance of their Fi parents. 
The differences can be illustrated as follows: The mortal- 
ity rates of the Koriyama Fi group and the Koriyama F 2 
group obtained from the Koriyama Fi adults were 53.0% 
and 79%, respectively, while the mortality rates of the 
Motomiya Fi group and the Motomiya F 2 group ob- 
tained from the Motomiya Fi adults were 31.2% and 
99%, respectively. As the mortality rates of the Okinawa 
Fi group and the Okinawa F 2 group obtained from the 
Okinawa Fi adults were 8% and 17%, respectively, the 
increase in mortality rates observed in the F 2 generation 
cannot be readily attributed to the transgenerational ac- 
cumulation of radiation effects. However, effects of the 
¥ l generation on the F 2 generation were discernable in 
the form of a small forewing size and its high correlation 
with the cumulative caesium dose, especially in females. 
It is important to note that a small forewing size can be 
induced not only through feeding of contaminated leaves, 
but also through long-term low-dose external irradiation 
in this species [7]. Furthermore, a small forewing size has 
been detected in field-caught butterflies from the polluted 
areas [7]. Therefore, the small forewing size observed in 
the present study implies that internal radiation exposure 
may have contributed to the high mortality rates recorded. 
Nevertheless, the Fi effect was small, and it could largely 
be overcome by the consumption of Okinawa leaves or 
amplified by the consumption of Koriyama or Motomiya 
leaves in the F 2 generation. That is, the effects were trans- 
generational, but largely non-genetic physiological ef- 
fects that could be overcome by dietary changes. 
Mechanistically, these outcomes may be explained by 
maternal effects, epigenetic effects, genomic instability, or 
similar phenomena. 

These results potentially conflict with our previous 
findings showing that morphological abnormalities, 
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Figure 6 Effects of the ingested caesium dose on the mortality rate in the F 2 generation, (a) A two-dimensional scatter plot. The ingested 
caesium dose in the F 2 generation (but not the F, generation) is provided on thex-axis. (b) A three-dimensional scatter plot. The ingested caesium 
doses in both the F, and F 2 generations are plotted. "F1" and "F2" indicate the ingested caesium dose in the F, and F 2 generations, respectively. 
"MR" indicates the mortality rate. Plots showing the same amount of F, ingestion are connected by lines. 



including wing colour pattern aberrations, found in the 
V 1 generation were heritable by the F 2 generation [7-9]. 
In this previous experiment, we used P-generation adults 
collected from the field in the spring of 2011. These P- 
generation individuals were directly exposed to all of the 
radioactive materials from the explosion of the NPP 
when they were larvae, which could have resulted in 
genetic damage. We then employed morphologically ab- 
normal Fj adults to obtain the F 2 generation using non- 
contaminated leaves. In the case of the possible genetic 
damage indicated in 2011, the larvae may have con- 
sumed radioactive materials adsorbed on the surface of 
leaves. The contributions of (3-ray-emitting nuclides in 
the earlier experiment cannot be ignored [14]. In 



contrast, the present study was performed in 2012, when 
the major radioactive species present was caesium, and 
the surface adsorption of radioactive materials was likely 
much less extreme. Therefore, the probability of genetic 
damage resulting from the ingestion of a contaminated 
diet was rather low in the present study. Both experimen- 
tal set-ups are reasonable for examining heritable effects 
that are likely based on either genetic or physiological (or 
epigenetic) changes. 

In the examination of forewing size, we noted a sex 
difference. Females showed a higher correlation with the 
ingested dose of radioactivity in both the ¥ 1 and F 2 gen- 
erations than males. The reasons for this sex difference 
are not clear. However, higher sensitivity of females to 
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Figure 7 Survival curves for the F 2 generation. The survival curves are segregated into two clusters of the experimental groups. For each 
curve, the localities of the F, and F 2 leaves are indicated, e.g., as Koriyama (F^-Okinawa (F 2 ). 



cold shock, a type of environmental stress, has been re- 
ported in this species of butterfly [4]. Hence, it is likely 
that females are generally more sensitive to stress than 
males of this species. 

The internal exposure of organisms living in a polluted 
area through a contaminated diet is unavoidable. Indeed, 
the possible biological effects of contaminated diets on 
humans after the Chernobyl accident have been docu- 
mented [15,16]. However, to the best of our knowledge, 
controlled experimental evaluations of internal exposure 
via a consumed diet are scarce in any organism in the 
field of radiation biology, with a few possible exceptions 
[17], although there may be some anecdotal or non- 
published studies that were performed following the 
Chernobyl accident. Thus, our experimental results con- 
stitute an important contribution to the research field. 
Accumulation of such scientific evidence is required to 
establish the low-dose biological effects of radiation [18]. 

Conclusions 

Our experimental system for detecting the biological ef- 
fects of internal exposure to radiation via contaminated 
diets revealed that low-dose ingestion likely causes death 
(i.e., high mortality and low survival rates), disease (i.e., high 
abnormality and low normality rates), and growth retard- 
ation (i.e., a small forewing size), at least in some indi- 
viduals of the pale grass blue butterfly. The effects of 
ingesting a contaminated diet on survival rates and fore- 
wing size were shown to be transgenerational. However, 
the effects of ingesting a contaminated diet in the F x 



generation could be largely overcome by ingesting a non- 
contaminated diet in the F 2 generation. It is likely that at 
least some of the biological effects of internal exposure 
to low doses of radiation are attributed to non-genetic 
physiological (or epigenetic) changes. 

Methods 

Ethics 

No specific permissions were required to collect the pale 
grass blue butterfly Z. maha and its host plant of O. cor- 
niculata in Japan. The pale grass blue butterfly is the 
most common butterfly in Japan, and its host plant is a 
weed that is often a target for eradication in gardening 
and agriculture. 

Butterfly rearing 

We followed standard rearing methods we have described 
previously [8], with a few minor modifications [4,7]. The 
egg collection procedure was also described elsewhere 
[7]. A pool of females was caught in the field in Okinawa 
(defined as the P generation) and used to collect eggs (de- 
fined as the Fi generation or the first offspring gener- 
ation). These eggs were divided into four (in the summer 
of 2012) or three (in the fall of 2012) feeding groups in the 
Fi generation. Therefore, the genetic differences among 
these feeding groups could be ignored. Eggs were similarly 
obtained from a different pool of females caught in 
the field in Okinawa and were then divided into three 
feeding groups. These groups were fed leaves from 
Okinawa, Koriyama, or Motomiya. These additional three 
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groups were employed in crosses with the first-generation 
groups to obtain the next generation (defined as the F 2 
generation or the second offspring generation). This 
process was necessary to avoid sibling crosses that could 
adversely affect the robustness of the offspring. To obtain 
the F 2 eggs, morphologically normal and healthy (non- 
abnormal) adults were selected from one Fi group and 



mated with healthy adults from a different F x group. In 
this paper, a group that consumed leaves from locality A 
in the Fi generation and leaves from locality B in the F 2 
generation is often indicated as the A-B group, but these 
butterflies all came from Okinawa lines genetically. All 
rearing experiments were performed in our laboratory, lo- 
cated on the Nishihara Campus of the University of the 
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Table 6 Pearson correlation coefficient r and p-values for 
the F 2 generation between the ingested caesium dose 
and forewing size 



Caesium dose 


r 


P 


Male, F 2 only 


-0.35 


0.44 


Male, F, + F 2 


-0.59 


0.17 


Female, F 2 only 


-0.58 


0.17 


Female, F, + F 2 


-0.93 


0.0028 



Ryukyus, Okinawa, except for the rearing of the Kanto- 
Tokai groups, which was carried out in our branch at 
Kichijouji, Musashino City, Tokyo, where the radiation 
level was maintained at 0.04-0.05 uSv/h, with the windows 
being closed at all times. 



Collection of host plant samples 

Bunches of the host plant leaves of O. corniculata were 
collected from various localities in the summer and fall 
of 2012, as shown in Table 1. The bunches were covered 
with wet paper, stored in a chilled container, and sent to 
the laboratory, where they were stored in a refrigerator 
until use. The leaves stored in this way were still fresh 
when they were given to larvae. The bunches were stored 
separately, depending on the collection locality, to prevent 
cross contamination. Feeding experiments using the leaves 
from the Kanto and Tokai localities (Kashiwa, Musashino, 
and Atami) and Okinawa were performed in the sum- 
mer of 2012, and feeding experiments using the leaves 
from the Tohoku localities (Koriyama and Motomiya) 
and Okinawa were performed in the fall of 2012. 




Figure 9 Morphological abnormalities in F 2 individuals. Morphologically abnormal parts are indicated by red arrows, (a) Dented left 
compound eye in a Koriyama-Koriyama sample, (b) Eclosion failure in a Koriyama-Koriyama sample, (c) Eclosion failure in a Motomiya-Okinawa 
sample, (d) Truncated right antenna in a Koriyama-Okinawa sample, (e) Eclosion failure in an Okinawa-Motomiya sample, (f) Truncated right 
antenna in an Okinawa-Motomiya sample. 
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Forewing size measurement 

We used an SKM-2000 digital microscope and the asso- 
ciated software SK Measure (Saito Kougaku, Yokohama, 
Japan) to measure forewing size in digital images. We 
measured the distance from the base to the marginal 
band at the end of the M 1 vein on the ventral side. Be- 
cause the feeding experiments involving leaves from the 
Kanto and Tohoku localities were performed separately, 
the Okinawa groups for the summer and fall experi- 
ments were both set at 1.00 in Figure 3, so that their re- 
sults could be shown in a single graph. 

Measurement and calculation of radioactivity 

The ground radiation dose at the collection sites was 
measured using an Aloka TCS-161 scintillation survey 
meter (Hitachi Aloka Medical, Tokyo, Japan) or an RAE 
Systems DoseRAE 2 dosimeter (San Jose, CA, USA) on 
the surface of the ground (0 cm high). Radioactivity 
measurements of caesium in the host plant and pupal 
samples were performed using a Canberra GCW-4023 
germanium radiation detector (Meriden, CT, USA) at 
the Instrumental Research Center of the University of 
the Ryukyus. Measurements were also performed using 
an Ortec GMX30 N-type HP germanium detector (Oak 
Ridge, TN, USA) at Nagasaki University to check the re- 
producibility of the measurements. The host plant samples 
were dried and burned to ashes (flameless combustion) be- 
fore being subjected to measurements, as described previ- 
ously [7,11]. Pupal samples were measured in a small 
columnar plastic container as indicated in a previous re- 
port [11], but the results were all below the detection limit. 

To calculate the amount of caesium ingested, we as- 
sumed that the larvae hatched on the sixth day following 
egg deposition, and they were subsequently fed Okinawa 
leaves for eight days, after which contaminated leaves 
from the different localities were provided. Thus, for the 
first 14 days after egg deposition, all larvae consumed 
Okinawa leaves. Based on previous data [7], we esti- 
mated that approximately 10% of the entire amount of 
leaves that a larva consumed came from Okinawa leaves, 
which was taken into account to calculate the ingested 
caesium dose. Furthermore, for simplicity, we assumed 
that the larvae ate all of the leaves required for subse- 
quent growth at once on the first day of ingestion of the 

137 134 

contaminated leaves and that Cs and Cs were re- 
leased at a 1:1 activity ratio on 15 March 2011 in a single 
burst from the Fukushima Dai-ichi NPP. 

Statistical analysis and graphics 

We used the statistical software R version 3.0.2 (R Foun- 
dation for Statistical Computing, Vienna, Austria), 
JSTAT version 13.0 (Yokohama, Japan), and JMP 11.0.0 
(2013) (SAS Institute, Cary, NC, USA). Linear regression 
equations were expressed with standard error and assessed 



with R 2 and j?-values. Data were also evaluated with 
Pearson correlation coefficient r. To assess forewing size 
differences, either Student's or Welch's £-test was per- 
formed, in accordance with F-test. Survival curves were 
examined by log rank test and Wilcoxon rank test. Data 
were graphically presented using Excel (2013) (Microsoft). 
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